Stable higher-charge discrete vortices in hexagonal optical lattices 
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We show that double-charge discrete optical vortices may be completely stable in hexagonal 
photonic lattices where single-charge vortices always exhibit dynamical instabilities. Even when 
unstable the double-charge vortices typically have a much weaker instability than the single-charge 
vortices, and thus their breakup occurs at longer propagation distances. 
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I. INTRODUCTION 

Some of the most spectacular experiments in the field 
of nonlinear light propagation in periodic photonic po- 
tentials relate to the properties of vortices and vortex 
flows [H . Self-trapped phase singularities of optical fields 
have been observed experimentally in the form of single- 
charge discrete optical vortices in square photonic lat- 
tices [3,[MI3l- addition, many of the theoretical and ex- 
perimental studies demonstrated that higher-charge dis- 
crete vortices are unstable, similar to the well-studied 
homogeneous nonlinear systems [l| . 

In this work, we study single- and double-charge dis- 
crete optical vortices in non-square periodic photonic lat- 
tices @, i, 0, S i, El- In particular, in the framework 
of the continuous nonlinear model of optically-induced 
lattices generated in saturable nonlinear media, we ana- 
lyze the existence, stability and dynamical properties of 
discrete optical vortices for the case of hexagonal optical 
lattices. We obtain the somewhat counter-intuitive result 
that double-charge discrete vortices in such lattices ap- 
pear to be far more robust and structurally stable than 
single-charge vortices. We verify this finding by demon- 
strating numerically the generation of a double-charge 
vortex with realistic experimental parameters. 

It is particularly important to highlight that although 
our results will be given with a view towards applica- 
tions in photorefractive crystals, they are not only rele- 
vant to that setting but also directly applicable to two- 
dimensional hexagonal waveguide arrays (e.g., in glass), 
showcased in recent experiments (see e.g., [ll| and refer- 
ences therein). Furthermore, they are likely to have di- 
rect implications to other areas of physics, such as Bose- 
Einstein condensates (BECs) in triangular lattices, the 
first experiments of which have just been realized [12l |. or 
even Debye crystals in dusty plasmas [31 . Another key 
aspect of the generality of our results is that they should 
also apply to honeycomb lattices. Hence, the findings 
presented herein have a bearing on two of the most fun- 
damental non-square lattice two-dimensional configura- 
tions. 
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FIG. 1: (Color online) Example of an always unstable single- 
charge discrete optical vortex for (3 = —0.76 (marked by a 
circle in Fig. [2} . (a) Intensity (top) and phase (bottom) ; (b) 
real (top) and imaginary (bottom) components; (c) absolute 
value of the corresponding Fourier transforms; (d) spectrum 
of the linearized equation displaying the linear instability of 
the configuration due to the presence of positive real parts in 
the eigenvalues A in the spectrum. 



II. THEORETICAL SETUP 

We study beam propagation through a self-focusing 
nonlinear medium in the presence of a two-dimensional 
hexagonal lattice by employing the continuum model 
with a saturable nonlinearity. To render our setting com- 
pletely amenable to the experimentally accessible regime, 
we use the theoretical model of a photorefractive nonlin- 
ear medium, which is known to exhibit strong saturable 
nonlinearity [8] . Polarization anisotropy of the nonlinear 
photorefractive response enables one to optically imprint 
various types of refractive index modulation (optical lat- 
tice) which can then be probed by an external beam [3] ■ 
Then the propagation of this beam in the presence of an 
optically-induced hexagonal refractive index pattern is 
governed by the following normalized evolution equation 



■ du ^ . 
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l + Ip{x,y) + \u\' 



= 0, 



(1) 



where u{x, y; z) is the normalized amplitude of the elec- 
tric field, z is the propagation coordinate, denotes 
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FIG. 2: (Color online) Family of single-charge vortices vs. 
propagation constant (3. Top: maximum real part of the lin- 
ear stability spectrum. Bottom: power P ~ J U^dxdy. 
The circle corresponds to the discrete vortex given in Fig. [1] 
The dashed line indicates another unstable branch which, for 
larger (3 bifurcates into different configurations. 
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FIG. 3: (Color online) Examples of (a) stable and (b) unsta- 
ble double-charge discrete optical vortices for /3 — —0.76 and 
/3 — —0.96 respectively (marked respectively by the circle and 
square in Fig. |3J . The layout of the panels is the same as in 

Fig.m 



the transverse Laplacian with respect to {x,y), D is 
the relevant diffraction coefficient, and 7 is the mate- 
rial parameter which is positive or negative depending 
on whether the nonlinearity is of focusing or defocus- 
ing character. The function Ip{x,y) = Ig\exp{ikx) + 
exp{—ikx/2 ~ iky^/3/2) +exp{—ikx/2 + ikyy/3/2)\'^ rep- 
resents the three-wave interference pattern that induces 
the hexagonal lattice. The lattice and beam intensities 
are normalised in units of the dark irradiance of the crys- 
tal. Throughout this work we use the following ex- 
perimentally realistic values for the system parameters: 
D — ZgX/ {AirnoxD = 18.015 (for laser wavelength in 
vacuum A = 532nm, and average refractive index of the 



FIG. 4: (Color online) Family of double-charge vortices vs. 
propagation constant /3. Top: maximum real part of linear 
stability spectrum [when nonzero, this denotes instability]. 
Bottom: power P — U^dxdy. The circle and square cor- 
respond to the stable and unstable discrete vortex configura- 
tions shown in Fig. [3] ((a) and (b) respectively). The dashed 
line indicates an unstable branch which, for larger /3 bifurcates 
into different configurations. 



2.36, Ig = 0.49, k = Air/Sd with 
and where the dimensions 
ys = l^m and Zs = 1mm 



medium no — 2.35), 7 
a lattice period d = 30/.tm, 
{x,y;z) are in units of Xs - 
respectively (see Ref. [7] for further details) 
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FIG. 5: (Color online) The top three panels (a) depict the 
evolution of a stable double-charge vortex configuration after 
a random perturbation with amplitude 5% of the initial am- 
plitude. The bottom three panels (b) are the evolution of a 
single-charge vortex configuration. In both cases (3 = —0.7. 
The color bar on the right provides a scale of the intensity 
(note the intensities of the single-charge vortex are lower rel- 
ative to (a) initially and saturated on this scale after break- 
up). 



We look for stationary solutions in the form 
u{x,y;z) — U{x,y) e'}q){il3z) expiiuKp), where U is real, 
f3 is the propagation constant, ip is the vortex phase, and 
m is the vortex charge. We solve the resulting nonlinear 
equation numerically, and the major results are summa- 
rized in what follows. 
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FIG. 6: (Color online) The same set of panels as in Fig. [T] 
except for a stable double-charge vortex in a honeycomb lat- 
tice. 
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FIG. 7: (Color online) (a) Input beam intensity profile rel- 
ative to the lattice (position of lattice intensity maxima are 
shown as rings). The intensity is given by the color-bar on 
the immediate right. Appearance of the beam sd, z — 60mm 
with different initial vortex phases (intensity not to scale of 
color-bar); (b) single-charge vortex, (c) double-charge vortex. 
Top panels: intensities; bottom panels: phase. 



III. NUMERICAL RESULTS 

We begin by considering the simplest six-site vortex 
structure, that of a single-charge (to — 1) discrete vortex. 
Figure [1] illustrates a typical example, while Fig. [5] shows 
the single-charge vortex linear stability (top) and power 
(bottom) as a function of the propagation constant in the 
semi-infinite band gap of the periodic potential. A posi- 
tive real part of an eigenvalue in the linear stability spec- 
trum leads to exponential growth of the corresponding 
linear excitation mode, and therefore to instability of the 
vortex. Somewhat surprisingly we can see in Fig. [DJtop) 
that the single-charge vortex has an eigenvalue with a 
positive real part across its entire region of existence, 
and therefore the single-charge vortex is always unstable. 

In contrast we find that double-charge vortices may be 
stable (see Fig.[21[a)), and even where unstable the insta- 
bility is weaker than the single-charge case (see Fig. S]) . 
In fact, as we can see in Fig. SJ^top) the double-charge 
vortex has a wide parametric interval where it is com- 
pletely stable (from —0.92 < (3 < —0.65), while outside 
this range it is unstable due to weak oscillatory instabil- 



ities (complex unstable eigenvalues, as evidenced by the 
spectrum in Fig. [3I^b)). We note that neither the single- 
nor double-charge discrete vortex families degenerate into 
a linear Bloch mode, as one can observe from the saddle- 
node bifurcation that occurs close to the edge of the first 
band of the linear spectrum in both Fig. [2] and Fig. |4l 
The various unstable single- and double-charge vortices 
which occur along the upper dashed branch in each figure 
are not discussed here. The typical evolution of the stable 
and unstable vortices is illustrated in Figs. [5l^a,b). Even 
though the single-charge vortex is lower in power than the 
double-charge vortex, break-up of the former into single- 
site fundamental discrete solitons occurs around z = 50, 
while the double-charge vortex has been propagated to 
z = 1000 with no sign of instability. 

To further our theoretical understanding, we employ 
a discrete model. In the latter the analytically tractable 
anti-continuum limit can be used, for which discrete vor- 
tex solutions can be explicitly constructed and a detailed 
stability analysis can be performed, as has been done for 
square lattices [l^. In such a setting we consider the 
six site configuration with topological charge to over the 
contour, which takes the form Uj = eKp(in(f)j) expiiPz), 
where (pj — 2t: jm/Q and j = 1, . . . , 6 for the six sites con- 
stituting the relevant contour. It is straightforward to see 
that this configuration yields non-trivial phase profiles 
for TO = 1 and to = 2. For these structures, according 
to the framework of [l5l | , the fundamental vortex will be 
unstable due to two double real eigenvalue pairs and a 
single real eigenvalue pair whereas the to = 2 configu- 
ration may be stable. These general results may also be 
physically understood as a consequence of the ID mod- 
ulational instability (MI) results )M\ along the ID (with 
periodic boundary conditions) six-site contour of the vor- 
tex. Such MI considerations predict that configurations 
where adjacent sites have less than a 7r/2 phase difference 
(i.e. a single-charge vortex) will be unstable, while those 
with more than a 7r/2 phase difference (the double-charge 
case) will be stable. 

It is important to point out here that, as the above dis- 
crete ID contour analysis suggests, our results can qual- 
itatively be extended to other cases where there exists 
a six-site closed contour, as e.g. in the so-called hon- 
eycomb lattice in which each index maximum has three 
neighboring maxima instead of six. A typical example of 
a stable double-charge vortex in a honeycomb lattice is 
presented in Fig. [6l Furthermore, extending our con- 
sideration of the ID six-site contour to the case of a 
defocusing nonlinearity, one can apply a so-called stag- 



gering transformation along the contour, Uj 



Substitution of this expression in the discrete equation 
transforms the model from defocusing to focusing (and 
vice versa). This amounts to translating the phase of 
every other node along the contour by tt and, hence, 
transforming an m = l(m = 2) focusing vortex to an 
TO = 2 (to =1) defocusing vortex respectively, suggesting 
a corresponding stability exchange. 
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IV. EXPERIMENTAL PROPOSAL 

Finally, we consider the generation of double-charge 
vortices and suggest parameters for their experimental 
observation. For our particular lattice parameters we find 
that generation of stable double-charge vortices is possi- 
ble over a wide range of input beam intensities and pro- 
files, at least within our isotropic medium approximation. 
We consider a Laguerre-Gaussian input beam with the 
profile shown in Fig. [3 kept as constant as possible as the 
input phase is changed, with maximum intensity ~ l.SIb- 
In the subsequent evolution we see break up of the beam 
into single-site discrete solitons if the initial phase corre- 
sponds to a single-charge vortex (Fig. [Tj^b)), while with 
an initial double-charge vortex phase we see stable gen- 
eration of the discrete double-charge vortex (Fig. EKc)). 
Output ai z — 60mm is shown, however we have seen no 
sign of instability in the generated double-charge vortex 
at a distance of z = 500mm. 

Based on the above considerations, we believe that in- 
puts of the type associated with m = 2 should be sus- 
tained during propagation not only by hexagonal crystals 
in photorefractive media, but also by two-dimensional 
hexagonal waveguide arrays (e.g., in glass), showcased 
in recent experiments jT]| . Importantly also, similar re- 
sults are theoretically expected and have been numeri- 
cally confirmed (data not shown here) to be valid in the 
case of honeycomb lattices in such media. 

V. CONCLUSIONS & OUTLOOK 

We have studied the existence, stability, dynamics and 
generation of single- and double-charge discrete optical 



vortices in two-dimensional hexagonal optical lattices in 
the framework of a continuum nonlinear model for pho- 
torefractive nonlinearity. We have found that, in con- 
trast to square lattices, double-charge vortices can be 
stable, while single-charge vortices are always unstable. 
Our main finding constitutes a general result for both 
hexagonal and honeycomb lattices that we expect to be 
verified experimentally. 

There are numerous directions along which it would be 
interesting to continue the present study. For example, it 
would be relevant to extend our analysis to the defocus- 
ing case for which our discrete theory predicts that the 
results should be inverted (i.e., that the m — 1 should 
be potentially stable, while the to = 2 will be unstable). 
Another direction of interest would be to attempt to gen- 
eralize such studies to genuinely three-dimensional, non- 
square lattice settings (e.g., in fee, bcc or hep crystals) 
and observe how realistic three-dimensional excitations 
may behave in these classes of models. Some of these 
directions are under present consideration and will be 
reported in future publications. 
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